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SUMMARY 
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This paper describes the use of linear rigid formats in COSMIC NASTRAN 
without DMAP procedures for the analysis of nonlinear propeller structures. 
Approaches for updating geometry and applying follower forces for incremental 
loading are demonstrated. Comparisons are made with COSMIC NASTRAN rigid 
formats and other independant finite element programs. Specifically, the com- 
parisons include results from the four approaches for updating the geometry 
using RIGID FORMAT 1, RIGID FORMATS 4 and 13, MARC and MSC/NASTRAN. It is 
shown that “user friendly" updating approaches {without DMAPS) can oe used to 
predict the large displacements and instability of these nonlinear struc- 
tures. These “user friend^," approaches can be easily implemented by the user 
and predict conservative results. 


INTRODUCTION 

The potential for high propulsive efficiency in the Mach 0.7 to 0.8 speed 
range has renewed interest in turboprop propulsion in recent years. Advanced 
turboprop concepts feature thin, highly swept and twisted propeller blades to 
achieve high efficiency and low noise at cruise. Deflections due to 
centrifugal loading, for blades of this type, exhibit a high degree of non- 
linearity requiring special analysis techniques not automatically available in 
COSMIC NASTRAN. 

As stated in the NASTRAN Theoretical Manual, "An important limitation of 
the automatic geometric nonlinear procedures provided with (COSMIC) NASTRAN is 
in the assumption that the applied loads from which the differential stiffness 
is derived remain fixed in magnitude and direction during deformation of the 
structure, and that their points of application move with the structure". 

This assumption is not valid for swept cantilever type structures lixe turbo- 
prop blades in a centrifugal force field where the direction ana magnitude of 
the force vector is determined by the displaced position of the structure. 

Highly swept turboprop blades in a centrifugal force field may experience 
an abrupt reversal of tip deflection or a large change in tip-chord angle as 
the rotational speed is increased. These sudden changes in position can 
severely affect aerodynamic performance and are referred to as instability 
points in this paper. The objective of this paper is to describe a study con- 
ducted on a representative advanced turboprop blade configuration. The 
selected configuration has a sweep angle of 60* which is at the high end of a 
range of sweep angles being considered for the advance turboprop concept. The 
purpose of the study was to investigate high tip displacements, inherent in 
these designs, using NASTRAN without DMAP alters. The study included the use 
of RIGID Formats l, a and 13 . Results obtained from these RIGID FORMATS are 
compared with each other and with other finite element programs. The math- 
ematical description of the various methods used is summarized in the appendix. 
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TURBOPROP GEOMETRY AND FINITE ELEMENT MODEL 

The model of the scaled turboprop blade shown In figure 1 is about 10 
inches long, has a tip chord of 2 inches, and a maximum chord at the hump of 
about 3.5 inches. The thickness varies from 1 inch at mid chord at the root 
to 0.040 inch at mid chord at the tip. The leading edge thickness varies from 
0.180 inch at the root to 0.022 inch at the tip. The trailing edge thickness 
varies from 0.077 inch at the root to 0.016 inch at the tip. The blade has a 
twist of 33.2*, a sweep angle of 60* and was assumed to be made of titanium. 
This blade was never built because of predicted high deflections and stresses 
at operating speeds. 

The NASTRAN finite element model is show in figure 2. This model consists 
of 423 grid points, 744 CTRIA2 elements and 2466 unrestrained degrees of free- 
dom. The model was restrained along the bottom of the shank. The middle node 
was restrained in all 6 degrees of freedom and the other six nodes were re- 
strained in 5 degrees of freedom with translation allowed in the direction of 
the spin axis. 


RIGID FORMAT 1 (STATIC SOLUTION) 

RIGID FORMAT 1 was used to investigate the turboprop geometric non- 
linearities in four ways: 

Case (1) (RFl-1) Updated Geometry/Angular Speed Difference - A centrifugal 
force field was applied to the model in increments of ISOO revolutions per 
minute. Nodal displacements and nodal forces (load vectors generated from the 
applied speed increment plus the summation of load vectors from previous 
steps) at each increment were stored in computer disk files. After each step 
grid locations were updated by adding the noaal oisplacements to the previous 
grid points and a new summation of force vectors was created by adding the 
applied load vectors to the previous summation of load vectors. The new set 
of grid locations defined the deformed model position for the next increment 
and the force vectors defined the load required to reach that position. The 
centrifugal force field applied at each step was equal to the square root of 
the difference of the new speed squared minus the last speed squared. The 
mathematical description for this case is summarized in section 1 of the ap- 
pendix. Seven NASTRAN runs were made from 1500 to 10 500 rpm. Tip rotations 
and positions were used to identify the onset of instability (fig. 3). For 
this case the large tip displacement between 4500 rpm and 6000 rpm (fig. 3(b)) 
was clearly indicative of the onset of structural instability. The change in 
tip-chord angle for this case (fig. 3fa)) was not as indicative of instability 
as the change in tip-chord position. The tip displacement also changes direc- 
tions as the angular speed increases from 4500 rpm to oOOO rpm (figure 3(b)). 
The tip displacement becomes progressively larger as the angular speed is 
increased to 9000 rpm and then changes direction again at 10 500 rpm indi- 
cating another region of instability. This method is straight forward and 
easy to apply using NASTRAN RIGID FORMAT 1. The predicted angular speeds 
causing instability are conservative (pessimistic) and, if accepted, will lead 
to design modifications to avoid the speeds causing instability. 

Case (2) (RFl-2) Updated Geometry/ Angular Speed - A centrifugal force 
field was applied to the model in increments of 1500 rpm. Nodal displacements 
at each increment were stored in computer disk files. After each step grid 
locations were updated by adding the nodal displacements to the previous grid 
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points. The new set of grid locations defined the deformed model position fc'* 
the next Increment. The tot?l angular speed was applied to this updated posi- 
tion. Thi> mathematical description of this case Is summarized In section 2 of 
the appendix. Seven NASTRAN runs were made from 1500 to 10 500 rpm. Again, 
tip-chord rotations and positions were used to Identify the onset of In- 
stability (fig. 4). For this case the onset of Instability can be noted from 
the large tip displacement and change In direction between 4500 rpm and 6000 
rpm (fig. 4(b)). Again, the change In tip-chord angle was not as Indicative 
of Instability as the change In position. The change In tip-chord position 
from 9000 rpm to 10 500 rpm Is more dramatic than In case (1). This method, 
also. Is straight forward and very ea ^ to apply using NASTP^ RIGID FORMAT 
1. The predicted points of Instability onset occur at lower rotor speeds than 
those using the method of case (1). Compared to the other two cases using 
RIGID F(MMAT 1, this case has the following advantages: (a) provides the most 
pronounced prediction of Instability, (b) requires only one computer disk file 
for storage of displacements, and (c) Is straight forttard and easiest to apply. 

Case (3) (RFl-3) Updated Ge^try/Angular Soe ed /Kestrainlng Forces - A 
centrifugal force fVelcl was applied to the model in increments of 1500 rpm. 
Nodal displacements and nodal forces (load vectors generated by the applied 
rotational speed) at each increment were stored in computer disk files. After 
each step, grid locations were updated by adding the nodal displacements to 
the previous grid points. The new set of grid locations defined the deformed 
model position for the next Increment. New force vectors were then created 
from the NASTRAN generated loads resulting from the current applied rotational 
speed. These force vectors were applied In the opposite sense at the next 
speed increment and represented blade Internal restraining forces (relative 
equilibrium with the deformed shape). The full rotational speed was applied 
at the next Increment. The mathematical description of this approach is sum> 
marized in section 3 of the appendix. Seven NASTRAN runs were made fran 1500 
to lO 500 rpm. Tip-chord rotations and positions were used to identify the 
onset of instability (fig. 5). As for the previous two cases, the onset of 
instability is noted by the large tip displacement between 6000 rpm and 7500 
rpm (fig. (5b)). The change in tip-chord angle was not as indicative of 
instability as the cnange in position. Tne point of instability is indicated 
by the change in tip-chord position from 6000 rpm to 7500 rpm (fig. 5(b)). 

This approach does not subject tne deformed turboprop to the full angular 
speed at the end of each increment. Therefore, the turboprop deforms the 
least compared to the previous two cases and, consequently, has tne smallest 
effect on the predictions for the vibration frequencies. Predictions of blade 
frequencies will be discussed later. 

Case (4) (RFl-4) Updated Geometry (Displacement Incrementl/Angular Speed 
Difference - A centrifugal force fieTd was applied to the model in increments 
of 1500 rpm. Nodal displacements and nodal forces (load vectors generated 
from the applied speed increment plus applied load vectors from previous step) 
at each increment were stored in disk files. After each step, grid locations 
were updated by adding the nodal displacement Increments (current minus 
previous) to the previous grid points and new force vectors were created. The 
new set of grid locations defined the deformed model position for the next 
increment and the force vectors defined the load required to reach that posi- 
tion. The centrifugal force field applied at each step was equal to the square 
root of the difference of the new speed squared minus the last speed squared. 
The mathematical description for this case is summarized in section 4 
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of the appendix. Seven NASTRAN runs were made from 1500 rpm to 10 bOO rpm. 
Tip-chord rotations and positions were used to identify the onset of in- 
stability (fig. 6). For this case the junp in chord position from i^OOO rpm to 
10 500 rpm (fig. 6(b)) is clearly indicative of the onset of structural in- 
stability. The change in tip-chord angle for this case (fig. 6(a)) is not at 
all indicative of instability. It shows possible angle-change leveling-off at 
higher rotor speeds. The tip-chord also changes direction as the angular 
speed increases from 3000 rpm to 4500 rpm (fig. 6(b)). The accompanying tip 
displacements become progressively larger as the angular speed is increased to 
9000 rpm and then changes direction again at 10 500 rpm indicating another 
region of instability. This method is also straight forward and easy to apply 
using NASTRAN RIGID FORMAT 1. The predicted angular speeds causing iiv 
stability are higher than case (1) or case (2) and about the same as case 
(3). This is the only case which predicts a possible leveling-off in tip- 
chord angle. However, intermediate points are needed to ascertain leveling- 
off in tip-chord angle with increasing rotor speed. , 

The nonlinear structural response predicted using RFl-4 is judged to be 
the most accurate of the various methods examined herein. The accuracy of the 
method relative to error estimates and relative to experiirental data remains 
to be determined. 


RIGID FORMAT 4 (RF4) (NONLINEAR SOLUTION) 

As stated in the NASTRAN User's Guide, "RIGID FORMAT 4 is a second order 
approximation to nonlinear effects of large deflections which computes the 
differential stiffness matrix by an iterative technique and treats the new 
matrix as a load correction. The internal loads are thus not linearly related 
to the applied loads". (See NASTRAN User's Guide). The mathematical descrip- 
tion for this case is sunmiarized in section 5 of the appendix. 

For each step the load was divided into six increments and the dif- 
ferential stiffness matrix was iterated until the value of tne ratio of energy 
error to total energy was less than 0.001. 

The iterative solution for the differential stiffness matrix indicated 
divergence for speeds above 7500 rpm. For the four load cases, 3000, 

4500, 6000, and 7500 rpm, the number of iterations >"equired to meet the con- 
vergence criteria was 2 , 3, 4, and 7, respectively. The results are plotted 
in figure 7. The points of instability are easily identified in figure 7(b) 
where the tip-chord displacement changes direction from 4500 to 6000 rpm. 

These results are analogous to those for case (1) with the exception that the 
tip-chord rotation increases less . th angular speed up to 7500 rpm. 

This method is automatically applied. Its shortcoming is that the new 
load increment remains fixed in magnitude and direction during deformation of 
the structure. However, judicious interpretation of the results can provide 
the instability information needed. 

RIGID FORMAT 13 (RFl3) 

(Normal Modes with Differential Stiffness) 

RIGID FORMAT 13 is a one step differential stiffness solution in which the 
applied loads are assumed to move with their points of application and remain 
fixed in magnitude and direction. The differential stiffness matrix is com- 
puted from the internal element loads resulting from the static solution. 

RIGID FORMAT 13 combines the differential stiffness with the original stif- 
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fness matrix to solve the eigenvalue problem and obtain the normal modes 
solution. The mathematical description for this case is summarized in section 
5 of the appendix. 

Results for tip-chord rotation and tip-chord position are shown in figure 
8. Though a change in tip-chord displacement occurs between 2800 and 8000 
rpm, it is difficult to Identify points of instability using this RIGID 
FORMAT. Note that the tip-chord displacement changes more slowly as the 
angular speed Increases above 5000 rpm. Compared to the methods described 
previously. RIGID FORMAT 13 predicts (1) tip-chord rotation which is somewhat 
more severe than RFl-4 and (2) tip-chord position change which is less severe. 
It appears that RF13 will yield reasonable solutions of the geometry non- 
linearities including frequencies with the possible exception of impending 
structural instabilities. 


COMPARISON OF RESULTS 

The tip-chord rotation is an Important design parameter. It affects tne 
angle of attack which influences the efficiencies, as well as noise and aero- 
elastic characteristics of the propeller. Predicted tip-chord rotations 
obtained as described previously are compared with each other and with those 
predicted by three other methods in figure 9. These three other methods are: 
U) ICP-Independent contractor program, section 6 of the appendix, (2) 
MSC/NASTRAN (McNEAL-Schwindler Company NASTRAN - new nonlinear geometry 
capability) obtained at Ames Research Center, private communication, and (3) 
MARC obtained at Lewis. It is noted that the MSC/NASTRAN and MARC predicted 
results were obtained using standard formats. These results are shown here 
for comparison purposes only and not to imply that this is tne best that can 
be done using MSC/NASTRAN and/or MARC. 

The tip-chord rotation curves predicted by Rfl-3, RF4 and MARC are almost 
identical up to angular speeds of 7500 rpm. Those predicted oy Rri-4, RF13 
and MSC/NASTRAN are about the same and somewhat higher. Those predicted by 
using the other methods, (RFl-1, RFl-2, and ICP) are, respectively, higher. 

The results of methods RFl-1 and RFl-2 follow almost identical paths for tip- 
rotation. However, the method of RFi-2 is much easier to apply. It should be 
noted that both tip-chord rotations and tip-chord positions are required to 
properly assess tne propeller structural dynamic behavior in view of possible 
instabilities at higher speeds. In addition, changes in tip-chord rotation 
affect directly the angle of incidence and, consequently, the efficiency of 
the turboprop. 


VIBRATION frequencies 

The natural frequencies of the blade in the design speed range are 
important design parameters which affect stability and life prediction. The 
results of RIGID FORMAT 13 (Normal Modes with Differential Stiffness) applied 
to the initial grid positions and to the updated grid positions of RFi-2 is 
shown in figure 10. These curves indicate that the instability point near 
9000 rpm and the resulting decrease in natural frequencies would be overlooked 
if position updating procedures are not applied to RIGID FORMAT 13 analyses. 
Also, the margins for interference with cyclic disturbances could be mis- 
interpreted and, therefore, adversly influence the turboprop design. 

Similar results comparing updated-position frequencies with path- 
independent frequencies are shown in figure 11. Tne path-independent frequen- 
cies decreases rapidly and at lower rotor speeds leading to overly pessimistic 
frequency predictions and, perhaps, to major design modifications. 
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CONCLUSIONS 


COSMIC NASTRAN RIGID FORMAT 1 can be used, without recourse to UMAP 
alters, to predict the nonlinear structural response of swept twisted pro- 
pellers. RIGID FORMAT 1 with geometry updating using total displacements at 
each step (RFl-1, 2) can be used to predict the onset of structural insta^ 
bill ties versus rotor speeds. These predictions are comparable to those 
predicted using methods which are independent of load path and over-estimate 
the nonlinear structural response. The path- independent method predicts 
frequencies which are overly pessimistic. RIGID FORMAT i with geometry 
updating using displacement increments (RF 1-4) can be used to predict the 
overall nonlinear structural response. This approach appears to be the most 
accurate. RFl-4 predicted results for tip-chord rotation are somewhat below 
those predicted using Rrl3 and those predicted by the McNeal-Schwindler 
Company NASTRAN. RIGID FORMAT 1 with geometry updating using total dis- 
placements and negative force (RFl-3) appears to predict nonlinear structural 
response comparable to those predicted by RIGID FORMAT 4 and MARC. Results 
predicted by RFl-3 are below those predicted by RFl-4 and underestimate the 
nonlinear structural response. 

00 ^ & 


0 



APPENDIX 


ORIGINAL PAGE IS 
OF POOR QUALITY 


MATHEMATICAL DESCRIPTION 

The equations describing the matheinatlcs of the several approaches used to 
obtain nonlinear solutions are summarized In Sections 1 to 6 below. The 
equations are presented In matrix form. This Is a concise description of the 
many degrees of freedom (DOF) that these equations represent. The notation of 
the variables used In the equations Is as follows: R denotes spatial 
position (geometry); u denotes deformation from an R spatial position; K 
denotes stiffness matrix; M denotes mass matrix; fi denotes rotor speed; 
a denotes Increment. The subscript ”o” denotes Initial position; the sub> 
script F denotes final position; the subscript 1 denotes the current 
Increment; the subscript j denotes a sumnation Index with limits defined in 
the given equation. All the variables are referred to the global coordinate 
system (x, y, z). The equations are sunnarlzed In the sequence used. The 
notation Is also summarized In section 7 for convenience. 

1. RIGID FORMAT 1, Case (1) (RFl-1) Updated Geometry/ Angular Speed 
Difference - The Incremental equations at the 1th step for this case are: 



At ao^ -s aop, {R^f 3 |Rp} (equations (b) and (1)) which Is the 
final deformed spatial position. The total deformation for this 
given by 


( 1 ) 

( 2 ) 

(3) 

(4) 

(5) 


case is 


jUpf - {Rp} - {Rq} 


The nonlinear geometric effects accumulate progressively in equations (1) to 
(5) with successive rotor speed increments aOi. Note tnat the geometry 
Is updated by adding the displacement obtained from equation (b). This ac- 
centuates the nonlinear geometric effects and leads to more severe nonlinear 
response and, therefore, conservative nonlinear structural response. 

2. RIGID FORMAT 1, Case (2) (RFl-2) Updated Geometry/ Angular Speed - 
Tna Incremental equations at the ith step for this case are: 


fi-l}- {«i- 2 } ^ {“i-l} 

(7) 


(S) 


/ 
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“t [s-i]'' h-i] {“i-i} 


( 9 ) 

( 10 ) 


The final deformed position |Rp} Is obtained from equations (10) and (7) 
at fi.| - Qp. The total deformation Is obtained from lRp| - |Rq| . The non- 
linear geometric effects accumulate faster for this case than l^r RIGID 
FORMAT 1, Case ( 1 ) as can be readily observed by comparing equations (5) and 
(10). RIGID FORMAT 1, Case (2) will predict more conservative nonlinear 
structural response than RIGID FORMAT 1 Case (1). However, RIGID FORMAT 1, 
Case (2) Is much easier to Implement and will Indicate possible unstable 
spatial positions at lower rotor speeds. 

3. RIGID FORMAT 1, Case (3) (RFl-3) Updated Geometry /Angular 
Speed/Restraining Forces - The Incremental equations at the 1th step for this 
case are: 




(11) 

“i " l“i-l * 


(12) 

[S-l] - [^i-1 


(13) 

[«l-l] - [ «i-l 


(14) 

I-.1 ■ [>,-.] ■■ fi 1 

“U’t 

[}] (IS) 


The final deformed position is obtained from eouatlons (15) and (11) at 
flj a op. The total deformation M rfI - Rof. This case 
differs from the previous two cases in the way the current displacement 
Is calculated in equation (15). The full centrifugal force field is not 
applied at the ith step. The first term in the braces in equation (15) is the 
full rotor speed at that step (o) while the second, which Is subtracted from 
the first, represents the cumulative centrifugal force (oj_i). This 
approach will (1) underpredict the geometric nonlinear effects on structural 
response, (2) predict higher rotational speeds for impending Instabilities 
and, as a consequence, will be nonconservative. Implementation of this case 
is comparable to RFl-1. 

4. RIGID FORMAT 1 Case (4) (RFl-4) Updated Geometry (Displacement 
Increment) /Angular Speed Difference. The governing equations at the ith 
angular speed increment (using the notation defined previously) are given by: 


- {«o ^ § ‘“J 


(16) 

(17) 
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( 18 ) 

(19) 



At LQ^ • AQp, |Up}ls obtained from equation (20). Subsequently 
the total deformation Is obtained from equation (21) and the final deformed 
spatial position is obtained from equation (16). Note In equations (16) to 
(20) that the nonlinear geometric effects enter the equations progressively 
with each successive rotor speed Increment a&i. Note, also, that the 
total centrifugal force field Is used at each R^_i position (equation (20)) 
Increment plus cumulative, as compared to RIGID FORMAT 1, Case (3) (equation 
(15)) where only the Incremental centrifugal field Is applied. 

5. RIGID FORMATS 4 (RF4) and 13 (RF13) Standard NA5TRAN FORMATS - The 
differential stiffness Is used In the solution for both of these formats. 

Using the previous notation, the governing equation for the differential 
stiffness Is (See the NASTRAN Theoretical Manual) 

hf - ['<0 ^ '^d l“o>] {‘‘K"o]{''o} * [*‘d (“o- Vl)] 

where K<j is the differential stiffness matrix and |Uo| is the initial 
incremental displacement used to form . It can be seen in equation 
(22) that the stiffness matrix [kJ » L Kq + Kh (Uq)j Is formed using 
initial position variables and that Op is applied to the Initial position. 

As a result, the intermediate spatial positions and, therefore, the 
progressive nonlinear geometric effects are not represented In either [k 1 or 
Q. Consequently, both RIGID FORMAT 4 and RIGID FORMAT 13 will underpredict 
the nonlinear geometric effects on structural response. These formats will 
predict higher rotor speeds for impending structural instabilities. Solutions 
obtained by using either of these FORMATS will be more unconservative than 
those obtained by using RIGID FORMAT 1, Case (3). 

6. Nonlinear Geometric Solutions Independent of Load Rath - The total 
deformation for his approach is given by 

h-l • [[^o] - “o K]]’* {“r- 

Equation (23) will predict infinite |Up} 

[[■•] - Hi ■ « 


[".]{«,}} Bll 

if 

m) 
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Equation (24) has the form of an eigenvalue problem ami, therefore, Mill 
cause equation (z3) to predict Infinite displacements when 

oj . ^ (25) 

or when the rotor speed (q) equals any of the natural frequencies (») of the 
structure. Equation (23) will overpredict the nonlinear geometric effects on 
structural response as well as Intending structural Instabilities. Cotv> 
sequently, the results will be overly pessimistic. 

7. Summary of the Notation used In this Appendix 

[K] global stiffness matrix 

[M] global mass matrix 

(R) spatial position (structure geometry) vector 

jui global displacement vector 

x,y,z global coordinates 

A 1 ncrement 

fi rotor speed 

w natural frequency 

Subscripts 

F final 

d differential 

1 Increment (solution step) Index 

j summation index 

0 Initial 


lu 
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Figure 1. - Turt)oprop sUgc end propeiltr. 



Figure 2 . - Turboprop propeller finite element model 
(4?J grid points, 744 elementsi. 
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(a) TIp-chort rotation. (W TIp-chortf poiition. 

Figure 3. - Rigid format 1. caseOI- updated geometry and forcaitinguMrspaaddifferenC' rasulh. 
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la' I ip-chorfl rotation. (bt Tip-chord position, 

figure 4. - Rigid format 1. case (2l - updated geometry/angular speed resufts. 
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ANGULAR SPEED, rpm 
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PROPaiER ANGULAR SPEED, rpm CHORDWISE, in. 

la) Tip chord rotation. (bl Tip-chord position. 

Figure S - Rigid format 1, case t?) - updated oeometrvtingular speed/restraining force re- 
sults. 
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(a) Tip-chord rotation. 


(bl Tip-chord position. 


Fiqurpo, - Rigid (ormat 1, case (4i - updated geometry displacement increment and forces/ 
angular speed difference results, 
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Figure?. - Rigid formet 4 solution results. 
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Figures, - Rigid format 13 solution results. 
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Figure 9. - Compirisons of predicted tip-chord rotation using various 
methods. 
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Figure 10. - Campbell diagram showing comparisons of predicted nat- 
ural frequencies. 
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Figure U. ■ Cimpbell diagram showing compiiisons of predicted natu- 
ral frequencies. 



